#include <iostream>
#include <cmath>
#include <cassert>
#include<fstream>
#include<cstdlib>
#include"grid.h"

extern bool debug;

using namespace std;

//------------------------------------------------------------------------------
//------------------------------------------------------------------------------
void Grid::make_edges ()
{
}

//------------------------------------------------------------------------------
//------------------------------------------------------------------------------
void Grid::compute_ls_coef ()
{
}

//------------------------------------------------------------------------------
// Preprocess the grid
//------------------------------------------------------------------------------
void Grid::preproc ()
{
   // compute outward face normal for boundary faces

   // compute outward normal for boundary vertices

   // create edges
   make_edges ();

   // compute outward normal for boundary edges

   // least squares coefficients
   compute_ls_coef ();

   // check least squares coefficients
}
